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Summary. This text is meant as an introduction to the theoretical physics of photon emission 
in hot and dense strongly interacting matter, the principal application being relativistic nuclear 
collisions. We shall cover some of the results and techniques appropriate for studies at SPS, 
RHIC, and LHC energies. 



1 Introduction 

The formation, observation, and study of the quark-gluon plasma (QGP) is one of the im- 
portant goals of contemporary nuclear physics. Put more broadly, this program pursues an 
ambitious agenda which consists of mapping out the properties of a relativistic many-body 
system possibly out of equilibrium. Another way to state this goal is perhaps more intuitive: 
One seeks to explore the phase diagram of QCD (QuantumChromoDynamics: The theory of 
the strong interaction). This motivation has driven an energetic experimental program for sev- 
eral decades, at facilities which have included the AGS (Alternating Gradient Synchrotron, 
at Brookhaven National Laboratory), and the SPS (Super Proton Synchrotron, at CERN in 
Geneva, Switzerland). This research is in full bloom today, with experiments being carried on 
at RHIC (the Relativistic Heavy Ion Collider, at Brookhaven National Laboratory), and soon 
to be started at the LHC (the Large Hadron Collider, at CERN). In parallel, lower energy ex- 
periments continue at the GSI (Helmholtzzentrum fur Schwerionenforschung, in Darmstadt, 
Germany) and will also be part of the FAIR facility at the same site. 

Many variables have been proposed as probes of the QGP 1 1 1, but we will concentrate here 
on electromagnetic observables, more specifically on real photons. This text will concentrate 
on the theoretical techniques to calculate electromagnetic emissivities: The main experimental 
results and procedures are discussed elsewhere t2l . 

Three decades of relativistic heavy ion physics have made abundantly clear that progress 
in this field is usually associated with the measurement and interpretation of more than one 
experimental observable: The many-body nature of this problem is manifest. Moreover, com- 
plementary observables are not restricted to those of a hadronic nature. In fact, as they interact 
only electromagnetically, real and virtual photons have the potential to probe the entire dy- 
namical history of heavy ion collisions. The mean free path of photons being much larger 
than typical nuclear scales, they will leave the interaction zone unscathed and therefore con- 
stitute penetrating probes of the systems under scrutiny. However, the absence of final state 
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interactions is both an advantage and an inconvenient, the latter being small production rates 
(a -C a s ). Furthermore, photons will be emitted throughout the entire history of the nuclear 
collisions: A quantitative interpretation of the asymptotic spectrum requires an equivalent un- 
derstanding of the signal, together with a thorough understanding of the background. 

This text is arranged as follows: A brief outline of the formalism to calculate electromag- 
netic emissivities is given, followed by a discussion of photon production at the SPS, then at 
RHIC. At this latter energy scale, we review the electromagnetic signature of jets, together 
with the information they reveal on the strongly interacting medium. It is important to insist 
here that this text should be considered an introduction to some of the principles germane 
to the theoretical evaluation of electromagnetic emission rates and yields; it is not a detailed 
review. 



2 Photon production: General principles 

A formalism that lends itself well to the calculation of thermal radiation of electromagnetic 
quanta is that of finite-temperature field theory. Quite generally, if one considers two hadronic 
states |f) and |/) that differ only by the absorption or the emission of a photon, the transition 
rate between those two may be written as (3) 

where tV is the proper four- volume, the S-matrix element is 5// = (f\ Jrf 4 x^ u (x)A' 1 (x)\i), Jfi 
is the hadronic electromagnetic current operator, and Ap (x) is the photon field. After writing 
a free vector field as A* 1 (x) = e' 1 (e, ,k ' x + e - '*'* ) / \/2(oV (e^ is a polarization vector and k = 
(o),k)), and noting that the matrix element is translation-invariant, one arrives at 

Rfi = ~Sv {2K)4 [ S (Pi + k -Pf) + S (Pi- k -Pf)} (f\Jn(0)\i)(i\M0)\f) (2) 

where one delta function is relevant for the absorption process, and the other for emission. 
Taking a thermal average over initial states and summing over final states, one may write the 
emission rate in terms of a spectral function: 

= (2^Mv(«,k) 

/- v (fi),k) = -^ e -^(2Kf8( Pi -p f -k)(f\M0)\i}(i\J v (0)\f) (3) 
z >,/ 

where Z = Y,i e ~^ K '' with K = H — [lN. Going from the current-current correlation function 
to the field-field correlator with the equation of motion d^d ii A v {x) = J v (x), and using the 
Schwinger-Dyson equation, one can write the emission rate as 

0)-=- =-7^-t~-5 Imn*((B,k) (4) 

cPk (2tt) 3 e/^-l MV 

where JI^, is the finite-temperature retarded improper self-energy 1 3 . 4 ] which, to leading or- 
der in the electromagnetic fine structure constant a, is equal to the finite-temperature retarded 
proper photon self-energy. Therefore 
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V* = - ~kf ^TT Imn Mv (». k) + ^(« 2 ) (5) 

From Eq. (O it is immediately clear why photon measurements in relativistic heavy ion col- 
lisions have the potential to reveal precious information about the strong interacting system. 
This equation is perturbative in the electromagnetic sector but is exact to all orders in the 
strong coupling. Photon measurements (and in fact electromagnetic observables in general) 
have the potential to carry direct information about the in-medium photon self-energy at fi- 
nite temperatures and densities: A quantity which is currently not calculable. Bear in mind 
however that this information needs to be convolved with the time evolution of the colliding 
system. 



2.1 Photons from thermal hadrons 



In a heavy ion collision at relativistic energies, a plethora of mesons are produced, from the 
light to the intermediate mass scale [5], To model their interaction, a Massive Yang-Mills 
approach can be used 1 6 1 to write down a Lagrangian: 



Jz? = ^TiDfiUD^W + \F£liM(U + W - 2) 



-4Tr 



-%Tr (d^UDvU^'F 1 ^ +D jl U^D v UF Rliv 



-AU Rli 



YTiF^ v UF Rflv U i 



(6) 



The fields are 



2i v 



= exp 



A*= 1(V M -A M ) 



F l'v = d^Ay R - d v A^ R - /go [a£' S ,Av R 



M 



m 2 K - 



(7) 



Note also that F x = 1 35 MeV, and that A,- is a Gell-Mann matrix. , and A^ are, respectively, 
the pseudo-scalar, vector and axial vector meson field matrices, respectively (not the A^ of the 
previous section). The unspecified parameters of the Lagrangian can be adjusted to reproduce 
adequate hadronic phenomenology. For example, the masses and widths of the p and of the 
a\ may be used as input (7). Importantly, the ratio of D— to 5— wave in the final state of 
the a\ — > pit decay is determined experimentally to be —0.062 ±0.020 |8|. Our choice of 
parameters yields a value of —0.099 [9J. 

A practical and intuitive alternative to the field-theoretic approach to the calculation of 
photon emission rates is provided by relativistic kinetic theory. For example, the differential 
emission rate, R, in the emission channel 1 + 2 — > 3 + y can be obtained either from the imag- 
inary part of the two-loop retarded photon self-energy at finite temperature, or equivalently 
from 
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(8) 



where .M is the invariant scattering matrix element derived from the Lagrangian introduced 
previously. The plus (minus) sign is associated with a boson (fermion) hadron final state. The 
distribution functions, f(E), are also chosen appropriately. 

A study of the interaction vertices will generate all possible tree-level diagrams for the 
generic reactions X + Y — > Z + y, p — > Y + Z + y, and K* — » F + Z + y. Here, the sets 
{X,Y,Z} represent all combinations of n,p,K,K* mesons which respect conservation of 
charge, isospin, strangeness, and G parity defined for nonstrange mesons. Using Eq. ([8j, the 
thermal production rates are obtained by coherently summing matrix elements in each channel; 
convenient parametrizations of the results are given in the appendix of Ref. |9|. Note finally 
that the use of effective Lagrangians still necessitates the introduction of hadronic form factors 
at moderate and high momentum transfers, as hadrons are composite entities. One expects a 
reduction in the 2-3 GeV region (in photon energy) that amounts to an approximate factor of 
4 t9l llOI . A useful and effective form of the hadronic form factor is a standard dipole term: 
F(t) = (2A 2 /(2A 2 — ?)) 2 , where t is the usual Mandelstam variable. For A, an appropriate 
value for phenomenology is that of a typical hadronic scale, A = 1 GeV 1111 . 

In addition to the contribution from a thermal ensemble of interacting mesons, such 
as the one described above, there also will be photon radiation from channels involving 
baryons. A starting point here may be the hadronic many-body calculations used for the 
interpretation of low-mass dilepton measurements 1131 . A key theoretical ingredient there 
is the in-medium spectral densities for the low-lying vector mesons, as can be understood 
from Eq. {5]l. The connection between the electromagnetic signal and the in-medium spec- 
tral density is made explicit through the famous Vector Dominance Model (VDM) [14|: 
Imn M v =Li=p ,m,(t> i m i I Si ) Int^Vv (z) . In effect, VDM identifies the hadronic electromagnetic 
current operator with the vector meson field operator. This approach has had considerable phe- 
nomenological success, and further refinements have also been considered 1151 . 

Another thermal source of radiation is introduced in the next section. After this is done, 
all purely thermal electromagnetic emissivities can be compared at a same temperature. 

2.2 Photons from thermal partons 

The calculation of the thermal emissivity of an ensemble of thermal partons proceeds simi- 
larly to cases discussed previously. The key ingredient is the retarded, in-medium, photon self- 
energy evaluated at finite-temperature. In this context, it is instructive to follow the chrono- 
logical order of the different theoretical investigations. Early calculations first proceeded with 
an expansion in order of the loop topologies. It is revealing to spend some time on the rate 
for photon emission obtained from one- and two-loop self-energies 1 10 1. This discussion fol- 
lows that in Ref. 1 10 1. Taking the imaginary part of the self-energies shown in Figure [TJ will 
produce the annihilation and Compton channels for photon production which involve gluons, 
quarks, and antiquarks. Note that if the photon is virtual those cuts will also yield processes 
like qq — > y* and &(g 2 ) corrections to it 1121 . all of which disappear when the photon goes 
on-shell. Taking the imaginary part of the QCD diagrams in Figure[T]yields the rate of Eq. ^ 
which, taking the Boltzmann limit for the distribution functions, may be integrated to yield 
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Fig. 1. Contributions to the photon self-energy in QCD, at one and two loops. 



,^ ,, Compton _ 

fR = 5 cm, 2 pa 
cPk 9 6k 2 



(9) 



We have considered u and d quarks, and an implicit assumption here is that CO = E\ + E 2 ^> T, 
where CO is the photon energy. The values of the two constants S\ and 82 are known 1 10], 
and depend on whether quantum statistics are used or their classical limit. Importantly, k c is 
an infrared cutoff which regulates the rates: As the exchanged quark goes on shell, the well- 
known logarithmic divergence manifests itself if k c — * 0. As is known from Hard Thermal 
Loops (HTL) analyses a propagator must be dressed if the momentum flowing through is 
soft on a scale set by the temperature T (3). Consequently the infrared contribution to the 
imaginary part of the diagrams in Figure[Tjwill come from the diagram on the left-hand-side 
of Figure[2] The blob on the quark line in the figure represents the HTL-resummed propagator, 
where thermal effects regularize the low momentum (infrared) behavior of the theory. For u 
and d quarks, this contribution may be written as 

n" v ( P ) = -6x \e 2 TZj (0Tr[<r(fc)^(p-*)y v ] do) 

where 



^«=W^ + ^«^, (ID 



the dressed propagator for a quark with momentum k, and 

n 2 

\ 

Ik 



777 

«£(*) = { -*D±|k| + - i ^[ft ) (WW) : Ffii(*o/|k|)] ) d 2 ) 



The sum over fcy in Eq. is a sum over Matsubara frequencies |3|, and the functions Qq 
and <2i are the Legendre functions of the second kind, namely 

Qi(z) = zQo(z)-l (13) 

The effective quark mass is given below. In the limit g — » the bare quark propagator is 
recovered. Note also that po will be analytically continued to CO. 
In this language, the bare propagator is written as 

^)=« + (^^-(^ d4) 



with 
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g±(k) ; 



-*0±|k|) 



(15) 



As we know from HTL, the cutoff is soft on a scale set by T: We may choose it to lie in the 
interval between m q and T, where m q is the thermal quark mass: m 2 = g 2 T 2 /6 (with N c = 3). 
Doing so, and adding the answer one gets to the ones from Eq. one arrives at 1 1011 161 



d 3 R 



5aa s 2 - 
9 2k 2 6 



" £ ln 



(0 



(16) 



The value of the numerical constant S3 is also known, and depends on whether Boltzmann or 
quantum statistics are invoked. It is remarkable that the final rate is cutoff-independent. 

Therefore, at this order in the loop expansion, a leading-order (in a s ) contribution obtained 
from the imaginary part of the diagrams on the left panel of Figure|2]and of FigureQ]emerges 
. Notice the effect of using the HTL propagator: If the thermal mass is set to zero, the rate 





Fig. 2. Contributions to the retarded, finite-temperature self-energy of the photon at one loop 
(left panel), and at two loops (right panel). The dark circle on the quark and gluon propagators 
represent the use of effective propagators, resumed using the Hard Thermal Loops effective 
theory. 



diverges. It is clear that HTL does its job of shielding the rate from the looming logarithmic 
singularity. 

Progress was made through subsequent evaluations involving higher order configurations 
in loop topology. In fact, taking cuts through the diagrams on the right panel of Figure [2] 
represented a somewhat involved calculation |17| that produced what first appeared to be a 
surprising result: For three colors and for 2 quark flavors, the two diagrams on the right-hand- 
side of Figure[2]contribute 



r m 32 5aa s 

Imii,, ^ = 

M 3tt 9 



7T hear 

a> 



(17) 



In other words, diagrams of a higher order in "naive power counting" (ff(a 2 )) (right panel 
of Figure |2), contribute parametrically at order a$ as do diagrams of a less elaborate loop 
topology (left panel, same figure). Technically, this is understood as follows: As the one-loop 
contribution produced a logarithmic singularity when the quark thermal mass {m q ) tended to 
zero, the contribution from two-loop diagrams also become singular, but linearly. Specifically, 
the factor which blows up is T 2 /m 2 . That, multiplied by a 2 (from the vertices), produces a 
net contribution which is 0{tt s ). Perhaps more physically, the photon formation time can be 
related to the virtuality of the quark emitting the photon in a tree-level representation of the 
emission rate. Owing mainly to an uncertainty principle argument, quarks with low virtualities 
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will take a long time to radiate photons. Long formation will then leave the process vulnerable 
to in-medium scattering, spoiling the eventual decoherence of the photon field. This is in 
effect a statement of the Landau-Pomeranchuk-Migdal effect (LPM) effect 1 18 1. Searching 
for all contributions at this order, a diagrammatic analysis revealed that, at leading order in 
a s , the only set of diagrams that contribute to the photon rate is obtained by summing over 
ladder diagrams 1 19], as that in Figure [3] The ladder contributions can be written as a linear 




Fig. 3. A ladder diagram contribution to the retarded, finite-temperature, in-medium, photon 
self-energy at order a s . The dark blobs represent effective HTL propagators. Diagrams with 
crossed rungs, for example, have been shown to be sub-leading 1191 . 



integral equation and re-summed, as shown schematically in Figured to yield a finite result. 
Unfortunately, the net rate has not been amenable to a closed analytical form, but a convenient 
parametrization has been given in Ref. |20|. 




Fig. 4. A schematic resummation of ladder diagrams. 



Now, the purely thermal sources of photons may be compared against each other, in order 
to get some feeling for their relative importance. There will be radiation from a thermal ensem- 
bles of partons, and that from a thermal ensemble of hadrons from the confined sector of QCD. 
This comparison is done for a temperature of T = 200 MeV, and shown in Figure|5] There, the 
rate for photon emission obtained from the hadronic many-body approach of Refs. Illll21|[22l 
by considering the square of the virtual photon invariant mass and taking the limit M 2 — > 0, 
is the solid line. The mesonic contributions discussed in Section [2~T1 including hadronic form 
factors, are represented by the dashed curve. The photon rate from thermal partons 1191 ob- 
tained from taking the imaginary part of Figure [3] is shown as the dashed-dotted line. At low 
momenta, the net spectrum is dominated by radiation from interacting hadrons. At a momen- 
tum scale of about 1 GeV, and above, the photons from thermal QCD shine on more brightly 
that those from competing sources. It is also interesting to note that the enhancement in rate 
provided by summing the infinite set of diagrams from Figure [3] give a net boost of a rough 
factor of 2-~2.5 over those obtained through Eq. J16t . At this stage, an optimistic summary is 
possible: If the energy/temperature is high enough, the photons from the plasma state of QCD 
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Fig. 5. A comparison of photon emission rates from different phases of the QCD phase dia- 
gram. The nature of the different sources is discussed in the text. From Ref. (9) 

will be a large component of the electromagnetic spectrum. Also, the early conclusion that the 
hadronic and partonic phases have very similar electromagnetic emissivities 1 10 1 might still 
be tenable: The QGP rates complete to leading order in a s are within a factor of two of the 
sum of hadronic contributions (solid and dashed curves). Bear in mind, however, that as the 
temperature and density evolve in time, a net yield of photons will come from an integration 
through different landscapes of the QCD phase diagram. Also, there will be other sources of 
photons other than the ones discussed up to now, some of which will have a thermal compo- 
nent, some of which will not. The final spectrum has the potential to contain all of them. We 
shall explore these other channels in sections to follow. 



3 Photons from nuclear collisions at CERN SPS energies 

The first measurements of photons at the CERN SPS were made by the WA80 (23) and CERES 
1241 collaborations. Unfortunately, only upper limits on an eventual thermal component could 
be set, owing in part to considerable systematic errors. Subsequently, the WA98 Collabora- 
tion identified a photon excess (over those produced by a superposition of nucleon-nucleon 
collisions) in central collisions of lead nuclei at a fixed-target energy of 158 AGeV [25 1. This 
data is shown in Figure [6] The collaboration compared its data to those obtained in different 
proton-induced reactions at a neighboring energy, scaled to that used by WA98. For transverse 
momenta such that pj > 1 .5 GeV, this data shows a clear excess beyond what is expected from 
proton-nucleus collisions. In the process of attempting to interpret theoretically the WA98 
data, one must keep in mind that one photon source that is unavoidable is that associated with 
the primordial (during the very first instants of the nuclear collision) interactions of nucleons. 
As we shall see later, considerable sophistication has been reached in the calculation of this 
contribution at high energies. However, at y/s =17.3 GeV, analyses are still subject to ambigu- 
ities. For example, the question of whether the partons colliding at these energies have some 
amount of intrinsic traverse momentum {kj) has remained open: Attempts to extract some 
meaningful values from experiments have remained inconclusive |27|. Then, the absence of 
a consistent and successful theoretical framework at the level of proton-proton collisions has 
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Fig. 6. The real photon spectrum spectrum as a function on transverse momentum, as measured 
by the WA98 collaboration. The pQCD estimates are from 1261 . 



made the later assessment of nuclear contribution such as shadowing and dynamical kj broad- 
ening (the Cronin effect) challenging. In this context, the WA98 collaboration has recently 
shown 1 28 1 an analysis of photons produced in proton-nucleus collisions at \[s= 17.4 GeV. 
However, this additional data has not the level of precision necessary to remove theoretical 
uncertainties. It is a fair summary to state that the WA98 data has generated a flurry of theoret- 
ical activity 1 29 1. We now concentrate on an approach which combines several of the elements 
discussed up to now |9|. 

The hadron and parton rates have been defined in earlier sections, and compared at a fixed 
temperature in Figure[5] As an accurate description of the direct photon in primordial nucleon- 
nucleon at SPS energies is still a matter unsettled, empirical scaling relations extracted from 
fits to experimental data 1301 are used. Starting from data acquired in pp collisions, one may 
naively extrapolate to a nucleus-nucleus (A + B) geometry with 

d 3 V rompt d 3 a/ p 

d 3 k <jf p 

where T^(b) is the nuclear overlap function, N co ]\ is the number of nucleon-nucleon colli- 
sions, and <Ji" p is the inelastic proton-proton collision cross section. The differential photon 
cross section is parametrized from experimental data [9|. Since intrinsic kj smearing effects 
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Fig. 7. Direct photon spectra from proton-nucleus collisions, scaled to central Pb + Pb colli- 
sions at SPS energies, with different values of the momentum broadening parameter (Akj), 
compared with experimental data [31 1. From Ref. [9|. 



should in principle be contained in this parametrization, the dynamical smearing is generated 
by fitting an additional component to the available proton-nucleus data. This is tantamount to 
folding the parametrized spectrum with a Gaussian distribution: 

/(**■) = -rf£#j tr * ,{A * ) (19) 

One then obtains proton-nucleus photon spectra, which may be scaled [25 ] to central Pb + Pb 
collisions at 158 AGeV; this is shown in Figure|7] We may thus put all the ingredients together, 
and convolve our rates with a fireball evolution model which is also able to reproduce low and 
intermediate invariant mass dilepton measurement results [32]. In this context, it is important 
to highlight the importance of the fireball's hadro-chemical composition in this analysis. At 
the SPS, a fit to several hadronic species has fleshed out a scenario where particles undergo 
chemical freeze-out before kinetic freeze-out |5|. In the temperature zone between the two 
freeze-out values, the conservation of particle number demands the increase of chemical po- 
tentials: These will play an important role on photon production. In the hadronic phase, for 
example, the reaction 7Tp — » Ttj is enhanced over the case with nil chemical potentials by a fac- 
tor e^^"' 7 . It is therefore important to get a good quantitative understanding of both chemical 
potentials and transverse expansion: The latter blueshifts the spectrum towards harder mo- 
menta. Bearing these empirical constraints in mind, one calculates the final spectrum shown 
in Figure[8] together with the WA98 experimental data. 

The ingredients used here consist of a (schematic) fireball evolution which is nevertheless 
consistent with the measured hadrochemistry and with the hydrodynamic expansion global 
characteristics, of a comprehensive set of photon emissivity from a gas of hadrons, of partonic 
emissivities complete to leading order in a s , and of an estimate of the nuclear Cronin effect 
deduced from measurements made in proton-nucleus collisions generalized to nucleus-nucleus 
events. It is fair to write that one of the main conclusions emerging from this study is that the 
direct photons measured at CERN SPS energies are consistent with hadronic properties from 



Photon Production in Hot and Dense Strongly Interacting Matter 



11 




12 3 4 



q, [GeV] 



10 




12 3 4 



q, [GeV] 



Fig. 8. Direct photon spectra from central Pb + Pb collisions at SPS energies. In the top panel, 
the contribution from different sources are shown separately: That from the hot ensemble 
of hadrons (long dashes), from the QCD plasma (dash-dot), and from primordial nucleon- 
nucleon collisions (short dashes). The sum is shown as a solid curve. In the bottom panel, the 
net signal is shown for different values of the momentum broadening parameter (Akj). The 
figure is from Ref. (9). 

he confined sector of QCD: The plasma contributes negligibly to the observed signal. The data 
does not demand initial temperatures beyond those found a moderate range: 200 < 7) < 240 
MeV. Even though there might be still some unavoidable uncertainties in theoretical analyses, 
this conclusion appears robust. 
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4 Photons from nuclear collisions at RHIC energies 
4.1 A theoretical interlude on jet quenching 

The advent of RHIC - with an order of magnitude higher centre-of-mass energy than what was 
accessible at the SPS - has presented the community with incomparable opportunities for the 
formation and the study of the quark-gluon plasma. This entire program has been an unquali- 
fied success, which continues to this day. Arguably one of the main discoveries at RHIC has 
been the measurement of the spectacular absorption of the high energy jets by the strongly 
interacting matter 1331 134-1 . In this context, it had been proposed early on that elastic final- 
state energy loss of fast partons would suppress jets in high-energy proton-proton events [35:]. 
Even though the medium created in nucleon-nucleon interactions turned out to be an ineffi- 
cient sink for hard partons, subsequent studies following this seminal work have suggested 
that medium-induced radiative energy loss could in fact dominate its elastic counterpart, and 
that "jet quenching" could be seen in nucleus-nucleus collisions |36|. A variable which char- 
acterizes this behavior and which uses the data in proton-proton collisions as a baseline is 
the nuclear modification factor, for a given nucleus-nucleus process (A+B) in a given class of 
centrality related to an impact parameter b: 

d 3 N AB /dyd 2 p T 

RAA{y ^ b) ^(N coll (b))d^/dy d 2 pT (20) 

where d 3 N AB /dyd 2 pj is the spectrum of a given species, say, measured in nucleus-nucleus 
events (A+B). The denominator of Eq. < |20b is the same spectrum measured in pp collisions, 
scaled by the number of inelastic collisions. Thus by design, if the nuclear process is only 
a superposition of what happens in proton-proton events, Raa = 1- The nuclear modification 
factor has been measured by different experimental collaborations, for many different hadron 
species, at RHIC. Results obtained by PHENIX for neutral pions are shown in Figure[9] These 
results are striking. At such values of transverse momentum, the yield is naively expected to 
be dominated by QCD jet fragmentation. If the jets have the same characteristics as they do 
in vacuum, then the nuclear enhancement factor would be unity. It clearly is not. In addition, 
all experiments involving nuclei performed prior to the RHIC era have shown enhancement 
not suppression: As mentioned previously this has been ascribed to the Cronin effect. This 
experimental fact has recently been studied in Ref. |38|, and the different theoretical models 
used in this context were discussed and compared in Ref. |39|. The current status of our 
understanding of the results shown in Figure|9]is that the exiting jets lose energy owing to their 
strong interaction with the surrounding medium created in the relativistic nuclear collisions. 
Using the framework of perturbative QCD, the energy loss by gluon radiation induced by 
interactions with the hot medium has been found to be the dominant mechanism |40). To 
quantify these statements, we continue to use the language of field theory at finite temperature 
and use the formalism of AMY (Arnold, Moore, and Yaffe), as applied to jet energy loss 

MM- 

In this framework, the hard gluon and hard quark (or antiquark) distributions, P g and P q , 
evolve in time as they traverse the medium. The evolutions are coupled, and proceed through 

MM 
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Fig. 9. The nuclear modification factor, Raa, for it , measured by the PHENIX Collaboration 
in Au + Au collisions at RHIC. The different panels (left to right, top to bottom) start by a 
minimum bias sample and proceed through increasing centrality classes or, alternatively, from 
central to peripheral collisions. The figure is from Ref. ]37|. 



where the integrals are over k and run from — °° to °°. The integration range with k < rep- 
resents absorption of thermal gluons from the QGP; the range with k > p represents anni- 
hilation against an anti-quark from the QGP, of energy (k—p). The quantities dT jdkdt are 
the transition rates for the various processes (q — > qg, g — * qq, g — > gg), and be ob- 
tained by taking the imaginary part of the parton self-energy in the soft thermal background. 
Technically, an integral equation can be written and solved 1421 1431 . In writing Eq. i2U , we 
used dl^g(p,k) = dr g g g (p,p—k) and similarly for g — » qq; the function in the loss term for 
S ~* SS prevents double counting of final states. The effect of the radiative energy loss on the 
parton distribution of an initially monoenergetic quark jet is shown on Fig. [10] The vertical 
lines show the average energies at different times. This plot shows how efficient the hot static 
(and here infinite) medium is in knocking apart the initial delta function. It also shows how 
rapidly the initial symmetrical profile will get distorted: It is important to work with the whole 
distribution, not just with the average energy. 

In addition to the energy loss (and possible gain) by the radiative emission of gluons (and 
also possibly absorbing them), collisional energy loss is also possible. This is often termed 
"elastic energy-loss": A slight misnomer as the total energy is of course conserved during an 
elastic collision, but we are concerned here only with the energy of a subset of the participants. 
A study of the qualitative modifications of the energy profile brought about by the inclusion of 
a collisional component is instructive. This exercise was done with an approximate treatment 
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Fig. 10. The evolution of a quark jet with initial energy £/ = 16 GeV propagating through a 
medium of temperature T = 400 MeV, where the vertical lines represent the values of mean 
energies related to the corresponding distributions. The figure is from 1481 . 



for the elastic energy loss, where only the drag and diffusion terms were kept |51|. We show 
here results with a more complete treatment of the collisional loss (and gain) 1 44 ] : This is 
done in Figure QT] It is important to reemphasize again that any observables that will depend 
on the details of the energy profile will in principle be sensitive to the details being shown in 
Figure [TT] a detailed analysis which includes those different aspects is being completed. It is 
however becoming clear that one should not neglect collisional energy loss. 

A complete study of the phenomenology associated with the jet energy loss scheme re- 
quires: 

• Modeling the initial-state geometry of the strongly interacting zone 

• Modeling the space-time evolution of this system 

• Modeling the initial distribution of hard partons 

Given these elements, the jets are then evolved using Eq. J2U against the changing back- 
ground provided by the evolution model. There are several ways the individual elements in 
the list above can be realized. An important constraint is provided by the need to reproduce 
as many complementary physical observables as possible; this again is certainly one of the 
lessons learned studying the physics of relativistic nuclear collisions. This statement actually 
highlights the connection between the physics of jet quenching and of photon emission: 

1. Hard quarks that lose energy to the medium and that will hadronize into jets will have a 
fragmentation component into photons and dileptons 

2. Jet quenching is a jet-medium interaction. This class of interactions also has the potential 
to generate photons, real and virtual 

Thus, hard jet quenching and photon production at moderate to high pj can not be treated 
separately any longer. Several of the approaches used to calculate the energy loss of hard jets 
to the strongly interacting medium can be used to illustrate the connection between photons 
and hadrons born out of jet fragmentation; a comparative of several of these formalisms was 
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Fig. 11. The evolution of a quark jet with initial energy Ej = 16 GeV propagating through 
a static medium of temperature T = 400 MeV, where the vertical lines represent the values 
of mean energies related to the corresponding distributions. The top panel shows the effect of 
including only energy loss (and gain) through elastic collisions. The middle panel shows the 
effect of including only the radiative energy loss and gain. The bottom panel includes both. 
The figure is from the model of Ref. 1 44 1 and is courtesy of Bjorn Schenke. 



recently done 1 45 1 in heavy ion collisions, where the time evolution model was the same in all 
cases. 

In AMY, an overall free parameter is a s = g 2 /4-K, where g is the strong interaction cou- 
pling constant. It is convenient to fit this to hadronic observables, like hadronic spectra or R/^. 
An example of such a fit is shown in Figure [12] where the energy loss has been kept of the 
radiative type, and where the collision zone undergoes 3D hydrodynamical evolution. Note 
that the experimental nuclear modification factor for neutral pions is correctly reproduced in 
central and midperipheral collisions with a s = 0.33. Naturally, keeping a s constant is an ap- 
proximation and this value should be regarded as an average effective coupling. Note that, in 
the calculation of Raa, the single 7t° distribution first needs to be correctly reproduced I47II48I . 
Note also that in a three-dimentional expanding medium, the transition rates in Eq. d — 1 b are 
first evaluated in the local rest frame of the fluid cell, then boosted to the frame where the cell 
has a flow parameter f3. Specifically, 



dr{p,k) 



dkdt 



(l-v ]er f3) dr{Rhk0) 

UKQutQ 



(22) 

local 



where Vj et is the jet velocity, ko and to are the local momentum and proper time, respectively. 

With some quantitative control over the energy loss of hard partons and over the time- 
evolution of their distributions, one may now go back to the evaluation of photon spectra. 
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Fig. 12. The nuclear modification factor for measured n s, shown for central collisions (upper 
panel), and for midperipheral collisions (lower panel). The evolution model is an ideal (non 
viscous) 3D hydrodynamics |46|. The figure is from Ref. 1481 , the data from Ref. |49|. 



4.2 Photons from Jet-Medium Interactions 



Recently, the in-medium conversion of a hard parton to a photon has been suggested as a new 
source of electromagnetic radiation |50|. Specifically, when a fast quark, antiquark or gluon 
traverses a hot QCD medium it may undergo an annihilation (q + q — » g + y) or Compton 
(<? + S ~~* 1 + 7) process with a thermal parton. The production rate of photons is then 



dR 



./■ 



e f \2 T 2 aa s 



8tt 2 



{fqiPr) + fq(py)) 



( AEyT 



21n - 



-Ca—Cc 



(23) 



where C\ and Cc are two numerical constants respectively related to the annihilation and 
Compton channel; Ca = 1 -916 and Cc = 0.416. As the jets will lose energy during their cross- 
ing of the hot and dense strongly interacting medium, it is vital to consider this energy loss in 
order to properly assess the importance of this contribution. 

As hard partons crossing the medium can emit gluons and lose energy, they may also pro- 
duce bremsstrahlung photons as they interact with the parton constituents. The bremsstrahlung 
rate is obtained by appropriately modifying that for gluon emission |43 1. The photon distribu- 
tion associated with the bremsstrahlung process is therefore given by 



dPy{p,t) 
dt 



dkP q q(p+k) 



dri^iv{p+k, P ) 

dkdt 



(24) 



Note that even though the photon energy can be high, its production rate is small. Photon 
emission can thus be calculated perturbatively, without energy loss of the fast parton to the 
actual photons. Convolving the photon distribution with the initial spatial distribution of jets 
will yield the final spectrum of bremsstrahlung photons. 
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4.3 Non-Thermal Sources 

A source of photons that has no thermal component is that of "prompt photons": Radiation 
associated with primary nucleon-nucleon interactions. In a heavy ion collision, these will 
contribute in the very first instants of the scattering process, they constitute an irreducible 
background to the photons associated with the thermal processes one hopes to study. There- 
fore, nucleon-nucleon collisions constitute an important reference and their understanding thus 
forms a baseline requirement for theory. 

The prompt photon production in a proton-proton collision is given by 

J 3 prompt V [,, , nw „A^d<J ab ^y X a X„ 

E —£ — - = L / dx a dx b g(x a ,Q)g(x b ,Q) < K y — — —=- 

V Xa^fs-pjey ) 

da ab^cd j ~\ 

The different K-f actors take into account NLO effects; they are evaluated here using the 
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Fig. 13. The spectrum of real photons measured in p — p collisions at RHIC. The two sets of 
data points correspond to two diferent analysis techniques. The low pj data (green squares) is 
from Ref. 1 65 1, the higher momentum data is from Ref. 1 66 1. 



numerical program from Aurenche etal. (63 1 . The photon fragmentation function, Dyi c is ob- 
tained from Ref. |64|. An important purpose of this part of the photon production calculation 
is not only a reliable background estimate, but also a check of the NLO photon formalism at 
RHIC energies. Setting the factorization scale Q equal to the fragmentation scale Qf and set- 
ting both equal to pj / v2 yields the photon spectrum show in Figure[T3l The high momentum 
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data correspond to direct photon measurements, while the softer photons are actually lepton 
pairs that are treated as internal conversion photons. Note that it is crucial to get an accurate 
representation of the photon yields in pp collisions, as this constitutes the denominator of the 
nuclear modification factor It is fair to write that the choice of scales yields a good agree- 
ment with the pp data, and therefore constitute an appropriate prerequisite calibration for the 
calculations to follow. 



5 Photon yields at RHIC 
5.1 ID Bjorken expansion 

Prior to using a more realistic scenario to be discussed below, in order to get some intuition it 
is useful to first benchmark the different photon sources using a simple time-evolution model 
of the nuclear collision. A ID Bjorken expansion |52 | is therefore considered. Also, the initial 
temperature is assumed to scale with the local density, producing a temperature evolution of a 
QGPas 



T(r,r) = T i 



TAV3 
T 



1/4 

(26) 



The jet evolves in the QGP until it reaches the surface or until the temperature reaches its 
critical value, taken here to be T c = 160 MeV, when a first-order phase transition occurs. The 
fraction of the QGP thereafter in the mixed phase is 

where rj is the ratio of the number of degrees of freedom in the two respective phases (partonic 
over hadronic), and xy is the time when the temperature reaches T c . Recall that in a Bjorken 
expansion, the initial time and temperature are related by 

T^--^— 1 dN (2M 

Ti% -mr Q ^i^ (28) 

where gg is the number of degrees of freedom in the QGP phase. The measured particle 
rapidity density, dN/dy, thus determines the lifetime of the QGP. 

A free parameter in the approach described so far is a s (or equivalently the strong coupling 
constant, g). This is determined by requiring the measured nuclear modification factor for 
hadrons to be correctly reproduced, as discussed previously. After setting the value of the 
strong coupling constant, the photon production calculation may proceed without additional 
parameters. Integrating the rates introduced earlier with the ID Bjorken model produces the 
yields displayed in Figure [141 

It is clear that identifying different sources of photons simply from the signal strength 
represents a challenging tasks, as the size of the background contributions has to be known 
precisely. An experimental observable with a different and unique qualitative behavior would 
be a genuine asset. In this context, the azimuthal asymmetry of the electromagnetic signal was 
recently proposed as an observable with such potential. More precisely, a sizeable azimuthal 
asymmetry characterized by a coefficient V2 can be expected for intermediate to largepr pho- 
tons produced in non-central high energy nuclear collisions. It is helpful to recall that the 
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Fig. 14. Photon yields at the energies of RHIC (left panel) and of the LHC (right panel). The 
time-evolution model used here is a ID Bjorken expansion. The difference sources shown are: 
Emission from a purely thermal source (th-th), photons emitted from a jet interacting with the 
thermal QGP (jet-th), bremsstrahlung photons (jet-bremss), and photons from a fragmenting 
jet that has exited the thermal environment. Note that the bremsstrahlung curves in Ref. 1431 
and those shown here are different [ 54 1 . 



anisotropy is defined in terms of the Fourier analysis of the particle yield dN/ pjdpjdtj), 
through the coefficients v n : 



dN dN 



pjdpjdfy iTZpjdpj 



1 +Y,^v n (p T )cos(n<j>) 



(29) 



The angle <j) is defined with respect to the reaction plane. At midrapidity, all the odd coeffi- 
cients vanish by symmetry and vi is left as the dominant one. Importantly, at the energies of 
the SPS and of RHIC, hadrons at low and intermediate pj always have a positive V2, owing 
to the hydrodynamic flow |55|. A good example is shown in Figure Q3] On the other hand, 
propagating jets should lose more energy when they are sent into a region where the medium 
is thicker, namely out of the reaction plane. As hadrons at intermediate and high pj get an im- 
portant yield contribution from jet fragmentation, a stronger jet quenching should lead to fewer 
hadrons with the appropriate value of their momentum emitted into this direction. Importantly 
this "optical V2" is not associated with flow but with absorption; the argument implies a pos- 
itive V2 for hadrons originating from jets. No quenching would mean a vanishing azimuthal 
anisotropy for hadrons born from a fragmenting hard parton. Measurements at RHIC, for var- 
ious hadron species, support this interpretation 1571 . Concentrating on intermediate energy 
direct photons, it is clear from previous discussions that photons from jets (i.e. from fragmen- 
tation and from jet-medium interaction) will compete with those originating from primary 
hard scatterings of partons in the nuclei. 

To further quantify these arguments, a numerical calculation is done for Au + Au colli- 
sions at RHIC (^/snn = 200 GeV). Using again a simple ID Bjorken expansion, the initial 
conditions are constrained by the measured particle densities. The overlap region of the two 
colliding contains the primordial hard partons which will hadronize after escaping. This initial 
almond-shaped profile is determined here using a simple hard sphere geometry. Three central- 
ity classes are considered: - 20%, 20 - 40%, and 40 - 60%; fixing the initial formation time 
at T; = 0.26 fm/c, the initial temperatures are respectively 7] = 370, 360, and 310 MeV. The 
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Fig. 15. The elliptic flow of unidentified charged particles in Au + Au collisions, at 130 A 
GeV, as a function of the multiplicity of the events (left panel). The vertical bars indicate 
the range of early hydrodynamical predictions. The right panel shows minimum bias data 
of the same nature, as a function of transverse momentum, pr- The different solid curves 
represent hydrodynamic calculations with different equations of state and different freeze-out 
temperatures. From Ref. |56|. 



results for the photon azimuthal anisotropy coefficient are shown in Figure \TE\ Indeed, pho- 
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Fig. 16. Photon x>2 as a function of transverse momentum, for three centrality classes in Au 
+ Au collisions at the top RHIC energy. The dotted line shows V2 for the photons originating 
from primary hard collisions and from the fragmentation of hard partons. The solid curve 
shows the anisotropy coefficient for photons from all sources. The dashed line shows the effect 
of considering all photons, but neglecting the energy loss of propagating jets. From 1 58 1. 
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tons from primordial hard scatterings and from the fragmentation of hard jets show a positive 
vaue of V2, owing to jet energy loss. Interestingly, the contribution from jet-plasma photons 
(in this transverse momentum window) is enough to make \>2 go negative. At lower values of 
pj, photon anisotropy should receive a contribution from hydrodynamic flow, as do hadrons. 
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Such a hydrodynamic calculation has been performed for low pj photons [ 59] and the results 
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Fig. 17. Photon as a function of transverse momentum, for Au + Au collisions at 200 A 
GeV, for a fixed impact parameter of b = 7 fm. Quark and hadronic matter (dashed line) con- 
tributions are shown separately, and the elliptic flow of 7T's and p's is shown for comparison. 
From Ref . g9) . 

are shown in Figure [T7] Clearly, the V2 coefficient for low pj photons has the potential to 
display a revealing transition from a behavior characteristic of hadronic matter, to one typical 
of partonic matter. A similar calculation exists for lepton pairs |60| and thus lies outside the 
topics to be addressed here. Note finally that a complete picture demands a calculation where 
the soft and hard photons are reconciled in a common framework, and treated consistently 
within a given realistic dynamical model. The results of such a calculation are presented in the 
next section. 

5.2 Photon yields and relativistic hydrodynamics 

Whenever computing electromagnetic observables in relativistic nuclear collisions, it is im- 
perative to use a model whose parameters have not been specifically adjusted to fit those very 
observables, but can address a wide class of hadronic phenomena as well. It is therefore im- 
portant to use "realistic" hydrodynamical simulations approaches, namely ones that reproduce 
many of the soft hadronic characteristics measured at RHIC. One such model is AZHYDRO. 
This ideal (i.e. where the shear and bulk viscosity coefficients are zero) hydrodynamic sim- 
ulation in 2+1D has been empirically successful at RHIC: It is able to reproduce the bulk of 
the data on soft hadron production for pj < 1.5-2 GeV/c. This is the case for hadronic mo- 
mentum spectra from central to semiperipheral Au + Au collisions, including the azimuthal 
anisotropies |56|. Apart from ideal hydrodynamics, AZHYDRO assumes a boost-invariant 
longitudinal flow velocity. For variables measured near midrapidity, this theoretical input is 
found to be appropriate. Going beyond this requires fully 3D hydrodynamic modeling. The 
different photon sources already introduced earlier are revisited briefly here, in the context of 
the hydro. 



22 Charles Gale 



For the photons resulting from jet-induced interactions, an important element is the initial 
momentum distribution of jets, at impact parameter b. This is obtained from 



dN ie ,(Q,b) _, r de i+ j^ k+ i et x;Xj 

! ' = T AB {b) £ / dx igi ( Xi ,Q) 8j (xj,Q)K— VVTT 

d l q T dy fr^J dt K(Xi - e>q T /^s) 



(30) 



The ^-factor includes and accounts for next-to-leading order effects. Here, a value of K = 1.7 
will be used, which is essentially transverse momentum independent 1611 . Importantly, as will 
be seen later, isospin effects are included at the parton level: 

8a(x a ,Q)= \^f a (x a ,Q) + ^-^U(x a ,Q)\R(x a ,Q) (31) 

where f a is the parton distribution function inside a proton. The second term, /„*, corresponds 
to the distribution function inside a neutron and is obtained by the appropriate substitution of 
parton species. Shadowing effects are included in the function R(x a ,Q) |62|. The factorization 
scale Q is taken as qj . 

The full predictive power of the relativistic hydrodynamic is now invoked, in connection 
with the evolution of the hard partons characteristics. The local transition rates in each fluid 
cell are taken from the solution of Eq. ( 121b . and these depend on the strong coupling con- 
stant a s , and on the local temperature T. This quantity is evolved in he longitudinally boost- 
invariant AZHYDRO with the EOS Q equation of state |56| which matches a non-interacting 
QGP above T c to a chemically equilibrated hadron resonance gas below T c at a critical temper- 
ature of T c = 164 MeV. A somewhat early fhermalization, Tj = 0.2 fm/c is assumed, to partially 
account for a pre-equilibrium stage 1591 . Further note that the transverse flow velocity at any 
vale of the space-time rapidity can be inferred from that at r/ = 0. Writing (3 = (/3j_,j3 z ), we 
have 

/3 ± (t,7],x ± ) = /3±(T ' T? = °' X±) , ft = (t,7 7 ,x ± )=tanhr ] (32) 
coshrj 

Integrating the rates with the time-evolution generated by the relativistic hydrodynamics 
produces observables one can then compare with experimental measurements. The hadronic 
variables produced by AZHYDRO are discussed elsewhere |56], and the photon spectra are 
shown in Figure[T8] The contributions from hard primordial scatterings are labeled "prompt" 
in the top panel and include photons from Compton and annihilation events, together with 
those from jet fragmentation. The radiation from the thermal components of the quark-guon 
plasma is shown ("Thermal QGP"), and so is that from a hot gas of hadrons ("HG"). Com- 
paring the sum of the different contributions (the solid curve) with the experimental data, it 
appears that the photons originating from jet-plasma interactions constitute a non-negligible 
portion of the signal in the momentum window 2 < pj < 4 GeV. The lower momentum data 
are dominated by thermal emission, while the radiation from hard collisions gradually takes 
over the whole spectrum, for growing values of transverse momentum. Most of the higher 
momentum data in the lower panel can be attributed to Compton and annihilation contribu- 
tions calculated from pQCD. The purely thermal contributions are subdominant there and the 
fragmentation contribution to the real spectrum is small, owing mainly to the energy lost by 
the propagating jet. Adding together all the sources produces a good agreement with the data, 
with no additional adjustment of the model parameters (hydro characteristics, a„, . . .). 

Another useful representation of the same data is the nuclear modification factor now 
defined for real photons as 
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Fig. 18. The spectrum of real photons measured in Au - Au collisions at RHIC. The top panel 
data is extracted following the same technique (identifying low mass dileptons with a virtual 
photon) as that used for the low momentum part of Figure [T3l and is for a centrality class of 
- 20%. The data set "PHENIX (1)" is from (65), while the data set "PHENIX (2)" is from 
1671 . The latter supersedes the former. The bottom panel is for a centrality class of - 10%; 
the higher momentum data there corresponds to a direct measurement and is from Ref. [68|. 
The different contributions are discussed in the main text. 
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we only consider y = in this work. Also, as advertised previously, the azimuthal anisotropy 
coefficient might help disentangle some of the photon sources. Both these projections of the 
data are examined. In what concerns R^, it is first useful to isolate some of the cold nuclear 
matter effects; this is done in the left panel of Figure [19] In these estimates, a considerable 
effect on the nuclear modification factor is caused by neglecting the jet-plasma photons. This 
amounts to a reduction of approximately 30% (at intermediate values of pj), as seen in the 
right panel of Fig. [19] The two extreme cases - where jet-plasma photons are present or not - 
bracket the experimental data; the current large error bars do not permit a choice. The apparent 
downward trend of the data is intriguing. Isospin contributes to this as noticed in Ref. |70|, and 
seen in the left panel. Notably, in the calculations presented here, the additional suppression in 
Rj^ originates from the fact that jets fragmenting into photons have lost energy. This consti- 




Fig. 19. The nuclear modification factor for direct photons produced in central Au + Au col- 
lisions at RHIC, calculated using 2 + ID hydro. In the pQCD contribution, the scale used is 
Q = \fPl : /2. The left panel contains no finite temperature effects and shows the importance 
of shadowing and and of accounting for isospin. In the right panel, the effect of a QGP and of 
varying the scale are considered. The dashed line shows results obtained without considering 
the photons originating from jet-medium interactions. From Ref. |69|. 



tutes an important conceptual link between hadronic physics and electromagnetic observables. 
The large values of the nuclear modification factor at low momenta pj < 6 GeV can be directly 
related to thermally-induced channels, and appear to exceed somewhat the central values of 
the experimental measurements: Note however that pQCD slightly underestimates the denom- 
inator in R^a- Correcting this would help bring the calculation closer to what is determined 
experimentally on this linear plot. Smaller error bars will also enable progress in identifying 
the respective sources, as well as eventual contributions from cold nuclear matter [71 1. 

Turning to the determination of photon azimuthal anisotropy, results of hydrodynamical 
simulations can also be compared with experimental data. Hydro calculations for low pj 
photons were shown in Fig. [T7] and Fig. [16] displays results for higher pj which however 
used a Bjorken model expansion. Both regions are unified using a realistic evolution scenario 
and a realistic initial profile. Considering 



Jq" d$ cos 2(j)d 3 N^(b)/d 2 p T dy 
J^d(j)d i Nr(b)/d 2 p T dy 



the results of the theoretical calculation are compared with the experimental data in Figure l20l 
One sees that the experimental anisotropy is small and in fact essentially zero for pj > 5 GeV. 
The v'2 coefficient associated with jet-plasma interactions is indeed negative, but its influence 
is greatly diminished as compared with earlier estimates of Figure[l6] associated with simpler 
dynamics and geometry. Future analyses with isolation cuts for the jets might however restore 
this sensitivity. 



6 Tagging jets with electromagnetic radiation 

It has become clear that differential variables will offer higher discriminating power than their 
one-body counterpart, in what concerns the dynamics of relativistic nuclear collisions. The 
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Fig. 20. The nuclear modification factor for direct photons produced in central Au + Au col- 
lisions at RHIC, calculated using 2 + ID hydro. In the pQCD contribution, the scale used is 
Q ~ \/P~t/2. The left panel contains no finite temperature effects and shows the importance 
of shadowing and and of accounting for isospin. In the right panel, the effect of a QGP and of 
varying the scale are considered. The dashed line shows results obtained without considering 
the photons originating from jet-medium interactions. From Ref. |69|. 



nuclear modification factor, Raa> f° r example, offers only moderate sensitivity to the under- 
lying probability distribution of energy loss |72|. In this context, correlation measurements 
could go a long way in constraining energy loss schemes. The yield of associated hadrons, 
given a trigger hadron, is sensitive to details of the evolution model for di-hadron correla- 
tions, much more so than R^a 1 73 74 1. Often termed the "golden channel", a very promising 
trigger of hard hadrons is a high pj photon |75|. The argument can be succinctly summa- 
rized as follows: At leading order, direct photons are exactly correlated to a recoiling parton, 
in a hard parton-parton collision such as one occurring in a Compton process or in an anni- 
hilation process. Modulo next-to-leading-order effects, measuring a hard photon on the near 
side yields valuable information about the far-side energy of the recoiling jet. Measuring the 
fragmentation hadrons then will directly contain info pointing to the lost energy. 

An observable amenable to the study of correlations, and which contains more informa- 
tion that one-body distributions is the so-called yield per trigger, a conditional probability 
distribution: 

P{Pt) 

In the above, P(pj) is the single photon distribution, is a hadron transverse momentum, 
and P(pj,pj) is the two-particle (photon-hadron) distribution 1761 : 

P(p r T ) = j^Jd 2 r^ AB (rx) 

x£ [dp J T P( F ! r )P(p r T \p J T ,r ± ,$), 
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Note that £?ab(?l) is the usual product of nuclear thickness functions, where the impact pa- 
rameter dependence is implicit. Explicitly, S^ab (b,r±_) = Ta (r± +b/2)7g (rx — b/2) / Tab [b) , 
and TAs{b) = J cfir^TA^ ^)Tb(t^ +b). The definition of the nuclear thickness function 
in terms of, say, a Woods-Saxon form is Ta(t^) = Jdz Pa(t_l,z). In addition, P{p J T ) and 

P{p T ,p T ) are the probability distributions for single-jet and back-to-back di-jet production. 
Studies of photon-hadron correlations make it useful to define the fragmentation function in 
terms of the yield per trigger: 

D M (zt,P T )=Pt P (Pt\Pt) 07) 

where zt = Pt/Pt- Finally, we define another nuclear modification factor in terms of this 
photon- triggered fragmentation function, normalized to that in pp collisions: 

, 7 , Daa(zt,P 7 t ) 
Iaa (zt ,P t ) = —— y-r (38) 

P>pp{ZT,P T ) 

We may now look at experimental data, with a theoretical approach which includes 

• All the photon sources discussed so far, but neglecting the purely thermal ones. We shall 
restrict our analyses to photon in the intermediate to hard sectors, where thermal emission 
can be neglected. 

• Both radiative and elastic energy loss (and gain) of the jets 

• The space-time evolution of the strongly interacting system, as modeled by relativistic 
hydrodynamics 

Before comparing with measurements, it is instructive to highlight the power of this theoretical 
tool by looking at variables that cannot be directly measured. For example, Figure [2T1 shows 
the probability distribution of away-side jets tagged by a 15 GeV photon on the near side, 
before the jet makes it out of the medium and fragments. Given a photon tag, this plot shows 
the jet probability distributions, for the cases where the jets are tagged by a fragmentation 
photon, by a photon which originated from an interaction of the jet with the hot plasma, 
and by a direct photon. This figure shows convincingly that the initial (i.e. before escape 
and fragmentation) jet has more energy on the average than the photon. Also, jets tagged by 
fragmentation photons dominate the high momentum part of the spectrum, whereas the jet- 
plasma photons are what tag the lower energy jets except immediately around the tag energy, 
where direct photons are dominant. Clearly, this information will eventually be encoded in the 
experimental observables, but the stage shown here is only accessible to the theory. Going a 
step further from Figure[2"Tl one may thus ask how this information will influence the measured 
hadron characteristics: This is shown in Figurel22l Given a 15 GeV near-side photon trigger in 
central Au-Au collisions at the RHIC top energy, the low pj yield per trigger for hadrons on 
the away side is dominated by jets whose tags are direct photons. In the higher pj domain, the 
fraction of hadrons coming from jets tagged by jet-plasma and fragmentation photons is large. 
Close to the trigger photon energy and higher, the fragmentation photons essentially represent 
all of the tags. 

It is exciting that data on photon-tagged hadrons are starting to appear, as the RHIC pro- 
gram presses forward. Figure [23] (left panel) shows the photon-triggered fragmentation func- 
tion , Daa{zt), as a function of the momentum fraction in mid-central Au-Au collisions at 




Fig. 21. Given a 15 GeV photon as a trigger in central Au-Au collisions, the initial jet dis- 
tributions are shown. The different contributions from fragmentation, jet-plasma, and direct 
photons (vertical line) are shown separately. From Ref. 1 77] . 




Fig. 22. Given a 15 GeV photon as a trigger in central Au-Au collisions, the yield per trigger 
of away-side hadrons is shown and is further separated in different contributions discussed in 
the text. From Ref. 1771 . 



RHIC together with the corresponding experimental measurements from PHENIX. The right 
panel of the same figure shows the ratio of that quantity to its pp equivalent. This ratio de- 
fines Iaa{zt)' Importantly, our approach predicts a flattening of Ijh for higher values of zt- 
Note that an accurate treatment of D pp is a necessary prelude to a successful interpretation of 
Iaa- This is the case here 1761 . The STAR Collaboration has also reported photon-triggered 
fragmentation function m easurements in central and peripheral Au-Au collisions |79|. These 
measurements, together with the appropriate calculations, appear on Figure [24] We find that 
the STAR measurements span the transition between the jet tagged by the direct photons at 
low zt, to those associated with jet-medium and fragmentation photons. Also, the STAR pe- 
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Fig. 23. (Left panel) The photon-triggered fragmentation function is shown as a function of 
zt, for mid-central Au-Au collisions, at RHIC. (Right panel) The photon-triggered Iaa as a 
function of the momentum fraction, for the same class of events. The data is from |78), the 
calculation from |76|. Several bins of hadron momenta are shown for comparison. 




Fig. 24. The photon-triggered fragmentation function is shown as a function of zt, for central 
and peripheral Au-Au collisions, at RHIC. The data is from |79|, the calculation from |76|. 



ripheral measurements are consistent with pp calculations. Clearly, correlation measurements 
such as the ones shown in this section have immense potential. On the theory side, a full NLO 
treatment will put the calculation on an even firmer footing, but results obtained to this point 
have been more than encouraging: Correlation measurements will represent a large part of the 
future. Beyond this topical introduction, it is useful to point out that the electroweak tagging 
of jets can be pursued with lepton pairs 1 80 1, and even with Z° bosons |81 1. Clearly, the full 
power of the techniques discussed here with those last observables will only be realized at the 
LHC. 
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7 Conclusions and outlook 

It is fair to write that the experimental program on electromagnetic probes in relativistic nu- 
clear collisions has been, and continues to be, a resounding success. Precision photon (and 
dilepton) spectra have been measured and analyzed at the SPS and at RHIC, and constitute an 
essential part of the program to be carried out at the LHC. In this context, predictions of pho- 
tons yields at the LHC have been made for the soft sector 1 9 1 , and also for harder photons [ 43 1 . 
Those two different kinematical regimes are summarized in Fig.[25]and[T4] These are but two 
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Fig. 25. Photon spectra predictions at LHC energies. (Left panel) The distribution for soft pho- 
tons is shown and is together with its different components: That from the quark-gluon plasma 
(QGP), that from a thermal hadron gas (HG), and that from primordial nucleon-nucleon colli- 
sions [ 9 1 . (Right panel) The panel shows the distribution of harder photons, with its primordial 
(NN), jet-medium, jet bremsstrahlung, jet fragmentation, and thermal components 1431 , see 
also El. 

predictions for LHC energies: A compendium of these can be found in Refs. 182 83 1. With 
its impressive increase in available energy, the LHC should venture deep into new territory. It 
is clear from the figure that the LHC's copious production of jets will feed directly into the 
photon channels. In fact, even the thermal processes will shine through, at the new collider. 
Finally, beyond the scope of this text is the story unfolding in hot systems with a high baryon 
density, such as those currently studied by the HADES dilepton experiment currently at the 
GSI, and later with the proposed CBM experiment at the same site. 

In what concerns photons produced in relativistic heavy ion collisions, the last decade has 
seen striking advances in both theory and experiments. We now have photon emissivities com- 
plete to order a s , for example. Electromagnetic emissivities have also been calculated for a 
clss of strongly coupled theories, utilizing string theory techniques [ 84 1. We have also reached 
a stage where photon measurements can help constrain the dynamical evolution models that 
are used in the modeling of relativistic heavy ion collisions. One of the reasons for this height- 
ened predictive power is the fact that both "rates" and "evolution models" have considerably 
improved in the sense of becoming more quantitative. Together with a sustained effort from 




in-med HG \ 

- QGP (T i= 845MeV) \ 
initial pQCD (pp) 
sum 
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the theory community, this progress would not have been possible in the absence of a superb 
set of experimental measurements of a large palette of complementary variables. Another fact 
that helps the electromagnetic phenomenology of heavy ion collisions is the fairly recent re- 
alization that hard hadronic probes and electromagnetic variables are no longer decoupled. 
No longer is it possible to treat jet quenching and photon measurements independently. Not 
only does this coupling provide empirical constraints for the theoretical treatments, but it also 
serves to enhance further the unity of the field. We have come a long way, and much remains 
to be done. 
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